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ABSTRACT 


This paper presents a package of nine programs for the 
TI-59 calculator. This package was developed as a solution 
BomewoO Problems. One problem involved expanding and modify- 
ing an existing set of programs; and a second problem involved 
developing five distribution approximating programs. The 
gay 20n to these problems represents a package with consider- 
able capability in computing confidence intervals, performing 
Meornesis tests and approximating distribution values. The 
@eoteribution approximations include inverse CDF values for 
the Normal, Chi-square, Student's t and F distributions, 
which allow the computation of confidence intervals without 
meane tables. 

The TI-59 proved to be a useful tool in solving these 
problems and demonstrated the capability of hand-held programma- 
ble calculators. The comprehensive set of user guides included 
in this programming package provides even the inexperienced 
user with a step-by-step introduction to this capability. 
Additionally, the methods used in preparing this programming 
package are directly eee Gerother calculators or 


computers. 
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PLL RODUCTLON 


Siem purpose OL thiS paper is to trace the development of 
a package of nine programs for use in the TI-59 calculator. 
Two of these programs compute confidence intervals for either 
the one-population or two-population situations. Another two 
programs perform hypothesis testing, again, for either one- 
Or two-population situations. The remaining five programs 
generate approximate distribution values for the Normal, 
Binomial/Multinomial, Chi-square, Student's t and F 
meetributions. 

The basis for this programming effort was a set of four 
T1-59 programs written by Professor P. W. Zehna for his 
personal use and later used in the classroom. Professor 
Zehna's programs also computed confidence intervals and 
performed hypothesis testing. However, his programs were 
not completely user-friendly, especially in terms of user 
guides; and they were dependent on obtaining some percentile 
meues trom standard tables. This paper then presents a 
Significant expansion in the scope of these early programs, 
The main thrust of this expansion includes a simplified and 
Standardized set of programs and user guides, while eliminating 
the dependence on distribution tables. 

This package of nine programs was designed for two types 
of users: the student, who might be asked to quickly solve 
Several very different problems in succession; and the working 


analyst, whose main concern is one specific problem which 





requires great accuracy. Of particular importance to both 
types of users is a detailed set of user guides (Appendix A) 
which includes sample problems. Many of these sample problems 
have been solved in two ways; the first way involves using 
program-generated percentiles; the other way requires input 

of tabled percentiles. While it is not absolutely necessary 
to look up tabled values when using the confidence interval 
programs, that option has been included, and should be used 
when increased accuracy is desired. Except for some cases 
involving small degrees of freedom, however, the accuracy of 
these first two programs 1S quite good as can be seen in the 
Sample problems of Appendix A. When tables are not available, 
the percentile values generated by the five distributional 
programs 1s outstanding (Appendix B) and can be used in 

lieu of tabled values. One very convenient feature of the 
distribution programs is the ability to provide values normally 
available only by interpolation in the standard tables. It 
Should be noted that the methods used to generate approxima- 
tions in the confidence interval program is slightly less 
accurate than those used in the distribution programs. This 
difference is due to the limited number of program steps 
available in the first two programs. Except for the two 
hypothesis testing programs, there is no requirement to use 
the applied statistics module in the TI-59. However, this 
module is required when performing hypothesis testing because 
of the large program size and the need to provide a signifi- 


Cance level with each test. Thus, it was possible to provide 





meeamred, yet simple, user guides to implement theoretically 
menrect and accurate programs, By providing the option of 
uSing either approximations or tabled values, both the student 


and the analyst can accurately solve a variety of problems. 
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Cee tURESOP THE SPROBLEM 


The programming effort presented in this paper was 
generated by two major problems. The foremost problem 
involved expanding and modifying an existing set of TI-59 
programs into a more user-friendly package. A second, 
closely related problem was adding the capability for generating 
accurate distribution approximations. The following discussion 
of the solutions to these two problems is a general overview of 
the particular solutions used. The specific methods and theory 


Seeoolution are left for later sections in this paper. 


A, USER-FRIENDLY PROGRAMMING 

User-friendly programming implies programming with the 
user's knowledge, ability, and hardware familiarity in mind. 
Im this light, a user-friendly program is a program which 
has significant capability, yet can be used by those with 
Sy a modest knowledge of either the calculator or the theory 
involved. Applying this definition as a framework for providing 
a user-friendly programming package resulted in four areas of 
effort. These areas are: standardized data entry, standardized 
solution procedures, maximized use of calculator capabilities, 
and improved user guides. 

i otandardized Data Entry 

The nine programs in this package all require some 

form of data entry. The data entry schemes for the confidence 


interval and hypothesis testing programs need careful 
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Standardization because of the similarity of the data and the 
Size of the data sets involved. There are three possible 
types of data in these first four programs. Data may be 
from a One-population sample, a two-population paired sample, 
Or a two-population independent sample. Any of these data 
can be easily entered in the form of summary statistics, 
if available; however, raw data requires some standardization 
between programs. The data entry schemes of the first four 
programs use Similar data entry subroutines which take 
advantage of the TI-59's data entry sequence. In each of 
these programs, the data entry subroutine is initiated by 
pressing [ |. The sequence of data entry then differs 
Slightly depending on the type of data being entered. All 
of these entry subroutines use a format which requires a 
to be pressed after each data point entry. This method 
saves one keystroke for each data point compared to the TI-59's 
two-keystroke, , entry method. Also, the jR/S 
key is very close to the numerical keyboard, compared to the 
2nd| key, thereby eliminating a source of data entry errors 
that frequently occur when using the Sequc weer 
The data entry schemes used for the five distribution 

programs presented no standardization problem. All of these 
programs, except the Multinomial, require only a few parameter 
entries. In the Multinomial program, which can accept as many 
as 35 pairs of parameters, the data entry problem was more 
Meeticuit. Each of the multinomial data points is stored 


Separately until computation begins. The data entry 
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sequence automatically repartitions the calculator to make 
room for this potentially large amount of data. 

Pitas tOnreca directly in registers has been somewhat 
Standardized between programs by uSing the same registers, 
where possible, for similar data. More information concerning 
the contents of data storage registers can be found in 
Appendix A, 

2. Standardized Solution Procedures 

A considerable effort was made to standardize the 
steps required in each problem solution sequence. in the two 
confidence interval programs (see Appendix A) problems are 
solved in three steps. In the first step, the calculator is 
repartitioned and the data are entered. The second step 
requires entry of either a percentage, in which case the 
program generates an approximate percentile, or a previously 
obtained percentile value. The last step involves the selec- 
meen and initiation of the proper solution subroutine. A 
Similar three-step sequence is used in the two hypothesis 
testing programs where the first step includes repartitioning 
and entry of test parameters. The second step involves data 
entry, and the third step, subroutine initiation. 

The five distribution programs contain only a limited 
amount of standardization due to the different nature of the 
programs. The Normal, Chi-square, Student's t, and F programs 
mae the a. a, labels for the same basic functions. 

The label is used in the Normal, Chi-square, and Student's 


(not F) to generate approximate density values. The label 


iS 





is used to generate CDF approximations in all four of these 
programs, where similarly, the label is used to generate 
inverse CDF approximations. The Binomial/Multinomial program 
Shares none of these standardized label uses. 

3. Maximized Use of Calculator Capabilities 

The ability of the TI-59 to repartition was the basis 
for the expansion in capability over that achieved by 
Professor Zehna's original programs. Repartitioning allows 
the addition of the program steps necessary for the inverse 
CDF approximating subroutines in the confidence interval 
programs. The added programming space was also used to com- 
pute one confidence interval estimate, for og? with LL known, 
which was not available in the original programs. Additionally, 
repartitioning makes possible an F distribution program and a 
Multinomial distribution program, 

However, repartitioning is not without its price. The 
larger programs now require three edges of the magnetic 
program cards, which presents a slight inconvenience in 
added loading and storage requirements. Other problems 
associated with repartitioning, such as inadvertent loss of 
program steps and unwanted or improper partitioning, have 
hopefully been eliminated from this programming package by 
extenSive validation with sample problems. 

4. Improved User Guides 

Preparation of improved user guides was a key element 

in making a user-friendly programming package. The complexity 


of the programs involved and the intended use by students 
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necessitated a departure from the standard TI-59 program 
record sheet. As can be seen in Appendix A, the improved 
user guides are organized with the user in mind. The general 
Outline for each user guide follows this pattern: 

fae Le roduction 

Eeeeeceneral Procedures 

c. Specific Procedures 

d. Additional Capabilities 

e. Labels Used 

f. Storage Register Contents 

g. Sample Problems 
The user guides for the five distrubiton programs have been 
combined to take advantage of the relative simplicity in 
using these programs. The answers to each sample problems 
are in a 10-digit format to provide a positive check on 
calculator output when working each problem. The confidence 
interval sample problems are solved in two ways to demonstrate 
the differences between using approximations and tabled 


values for the required percentiles. 


B. ACCURATE DISTRIBUTION APPROXIMATIONS 

The second major problem addressed in this section 
involves generating diStribution approximations for both the 
confidence interval programs and the distribution programs. 
The capability to compute accurate distribution approximations 
provides a new dimension to this programming package by 


eliminating the need for standard distribution tables when 
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solving confidence interval problems. A drawback of this 
new capability is the time required to obtain some values 
from the approximating programs. The time required by the 
confidence interval approximations is not nearly as long as 
Mmemeetor the distribution programs; but then, the quality of 
the approximations is not as good either (see Section III). 
Appendix B contains comparisons of tabled values with 
both the confidence interval program approximations (Type I) 
and the distribution program approximations (Type II). The 
actual probability values used in these comparisons were 
obtained using the distribution programs and can be regarded 
as being very close to the actual probability achieved. The 
missing values in the inverse F comparison are due to the 
inability of the Type I apprcximation to generate inverse F 
values when either of the degrees of freedom parameters is one. 
Also, while only selected approximations are listed in the 
comparison, the inverse CDF approximations are of nearly 
equal quality over the entire range of the appropriate 
function. The Type I approximations displayed in Appendix B 
are generally not as accurate as Type II approximations in 
terms of actual probability achieved. Only the inverse CDF 
meproximations for the Chi-square, t, and F distributions are 
presented in Appendix B since all the other approximations 
Which are available in the distribution programs duplicate 
table entries. These other approximations include probability 
Values, CDF values, and various other distribution values 


(see Appendix A). The superior quality of these approximations 


ho 





can be attributed to using the same approximating methods 


Used in the TI-59 Applied Statistics Module (see Section III). 


i 





ee elite ORY, 


This section on theory is presented as a background for 
the solution methods used in the accompanying programs. As 
such, it is not intended to be a primer in statistics. 
Instead, this section should be considered as an intermediate 
level derivation of the specific statistical methods used in 
Eeeoeramming. For a more basic explanation of this material, 
the references cited within each subject area, or equivalent 
texts, should be consulted. 

There are five subject areas discussed in this section. 
First, the theory used for estimating confidence intervals 
in the first two programs is discussed. Next, the hypothesis 
testing theory necessary for programs three and four is dis- 
cussed. And lastly, the methods used in all nine programs 


MemOotain approximations to distribution values are discussed. 


A. THEORY FOR ONE-POPULATION CONFIDENCE INTERVAL ESTIMATION 
The derivation of theoretical interval estimates will be 

done in the same order as these estimates appear in the User 

Guide for Program 1 (Appendix A). Most of these derivations 

use the pivotal-qualtity method to obtain confidence intervals 

(1,u) [Ref. 1: pp. 379-389]. Other methods used here will 

be discussed in slightly more detail, but will still be brief 

compared to the referenced texts. When forming a C.I. with 

a nonsymmetric distribution the interval will represent an 


equal tails solution, where equal tails implies 
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Meee |) = P(X > uj = (1-Y)/2. This method does not provide 
fmeeshortest C.I. for a given Bee however, this method is 
Bemmoniy used for its ease of computation [Ref. 1: p. 382]. 
Regardless of the method used, the resulting C.I. given in 
this section by (1l,u) represents the formula used to calculate 
interval estimates. 

mac. ft. for Normal Lk with oF Known [Ref. 2: pp. 77-80] 

Assuming that X is distributed NCU, C2) and using 


(K- [L)/(O/4f 0) as the pivotal-quantity, a 100°/'% Cees 


constructed thus: 


P ees -py/2 =< (X~- (0 )/( O/s[n) < (+ 7)/2| 2 of 


Substituting z = 21+ )/2 = eel ayy and simplifying we have: 


P Iz ~ 20/Afn <= <=: X + 2 O/T - 9 
(lu) = (x- 20/fn , x + 20//n) 


feo. lt. ror Normal Ul with oe Umimewmelkhet. 22 p. -o0- 


Neue Spe 277° Ref. 1: p. 381] 
Assuming X is distributed NCL, CTF ) and using 
weil )/(s/./n), Wiheners  disStrrputed ti(n-1), as the pivotal - 


quantity, a 100°¥% telecom constructed thus : 


P Pseyyyele) = (R- LL) /(s/yfn) < tasy/2| =. 


eapstituting t = ebay) cn? = aclaiegiayy ice 1) and simpli- 


fying we have: 


eS 





Px - ts/y~n= [lak + ts/yfn|= 7. on 
ie) = (x - ts/rfn , x + ts /sfn). 


Mme rl. tor Bernoulli p [Ref. 4: pp. 376-381] 
The pivotal-quantity method does not work for the 
Bernoulli case and another method will be briefly developed 


meer This method starts with two numbers a and b such that 


pla< t=] Fas 


From these limits we have: 
P {X = al =e) 2a P {X = »| , or equivalently, 
1 ba = nal = (1-Y)/2 = P pa > na 


Now )X, is distributed Binomial (n,p) which, for k <n, can 
be explicitly related to the incomplete Beta function and 


Memce to the F distribution to yield: 


(nX + LIP Vy /26%1 V2) 


Bo (a= RD CAR LIF yy 70070 ¥9) 
nk 
“(nk + (n= nk # WP pay y a (vgt2e¥-2) 
where v, = (2nX + 2), and v, = (2n - 2nX) 


1 and u form a 100Y % conservative random interval thus: 
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Pll <p <u =e 


The outcome (1,u) is a conservative C.I. in the sense that 
Memeeeiis discrete distribution the confidence that (1,u) 
contains p is at least 100YS%. 
meeec.l. for Normal o< TARtciKnOWmeP Ret. 5. penl27 5] 
Assuming that X is distributed NCL, OM ) and using 
Y((x,-)/)*, which is distributed X%(n), as the pivotal- 


quantity, a 100 Y% equal tails C.I. is constructed thus: 
ee 2 
PIV yyatn) = LOG MIT" = Xray ysl] = 


Substituting q,; = X12)» and q, = X*aeyys2™ 


and Simplifying we have: 
mie) = (ix. - njl“){1/q, . 1/a,). 


meee C.1l. tor Normal Ge with LL Omimowne (Ret. dim. 56 22 


Memes: Dis, 2e7 | 
Assuming that X is distributed NCU, O*) and uSing 
meijs-/o-. which is distributed X“(n-1), as the pivotal- 


quantity, a 100 Y% eGuidtmetai is C.i., 1S constructed thus: 


Z 2 Z Zz 
P x (1-y)seb2) = Cais en <z x eyy/2t2)] = aye 





Supotituting q, = X ay y/2-0, and ce Xap /2O-) 


and simplifying we have: 
Pp (n-1)s“/q < (0° < (n-1)s“/a,| = aa Gss 
Zz i f 


{1,u).= (n-1)s“(1/q, , 1/q, ). 


fee. 1. tor Exponential Our Pema 5 «oe Zoe 


Assuming that X, x XO are exponential random 


220° 
variables with parameter » and using 2 Ank, Tein eet 


tributed Ren). as the pivotal-quantity, a Oona equal 


meets) C.1. is constructed thus: 
; PX? ayy fol 2nd/2nk <2\rk < XM" (aeyy/2l2)] =). 


Simplifying we have: 
iP DP ry) /22n)/2nk ~ d < Xap yjel2n/2n| =, or 
(Lu) = (XE (ayy jol2nd/2nk + X° perv ysy(2n)/2nk). 


The C.I. for the mean time to failure ( Sie 1/)) is 


Senstiucted by inverting the above interval to yield: 


(Lou) = (2n%/ X74 my yyol2n) 4 2nk/X cy vy yol2n)) 
[Ref. We p. 382] : 


Dee 








Pee PHeEORY FOR TWO-POPULATION CONFIDENCE INTERVAL ESTIMATION 
The confidence interval estimates for two-population 
Situations are discussed here in the same order as they 
appear in the User Guide for Program 2 (Appendix A). All 
af these estimates use the pivotal-quantity method discussed 
earlier [Ref. 1: pp. 379-389]. As in the one-population 
case above, the C.I. given by (l,u) represents the formula 
Weed tO calculate the interval eStimate. 
ime c.L. For Bernoulli Py "Py tOmebanve ml and nm Retom 2; 
pee 249) 
Large m and n means NPy» mpy ,n(1-p,), and m(1-py) all 
Sreater than five. With this condition met and assuming that 


X and Y are independent and normally distributed, 


oe - Y) = (py - Dy ) 


V sy°/n + sy“/im 





Which is distributed approximately N(0,1), is used as the 


pivotal-quantity. A anon aa Clee hs eCOMstructed tnus ; 


Caan (p, - Py) 


P |-z 
EQ) 22 eS z 
r eee) a + ¥(1-Y)/m Lay aen i 
Substituting ec = X(1-X)/n + ¥(1-Y)Mm , and simplifying yields: 


Pp E-F- ce sy yy ery =< X-Y+t Hey i/ mye Ot: 
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ru) = (x--y = ely 7/2 ’ any a AeGy 2 


mee C.1, for Normal Hy- Uy LO cranam cr aired [Rein 2- 


peal 23 |] 


Assuming that X and Y are normally distributed and 
letting D = X - Y, and using (D- LL) /(Sp/ Yn), Hegel 3S 
distributed t(n-1), as the pivotal-quantity, a 100 3 Gaal S 


eometructed thus: 


Pas tituting iy = Hy - Uy, and t = Zar aye. 


Beg y)/2'""1) and simplifying yields: 


. 5 t55/ fh = Hy-My = B+ e5p/Yn] = 7 gO 


ya, *-( a, )/n 
mel,u) = (4- ts)/fn,d+ts,A/n), where Sp= ——~Taypyp 


and d. = 


a ey 
- ue 2 _ 2 
fees. Lt. ror Normal [Ly- [ly with Or ~ Cy Gh 


fRetae2: p. 123] 


Assuming that X and Y are independent and using 
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Meee) - (Lily - Ly) 


5 Mvilenwiowdiatribited t£(mtn-2} , 
(1/m + 1/n)S., 


as the pivotal-quantity. Here and elsewhere in this paper, 
n and m represent the number of data points in X and Y, 


Beepectively. A 100°Y% C.1. is constructed thus: 


(Z-¥) - ( {Ly- Uy) | 
ee )/2 0 a2 ary amr) |* 7- 
p 


Subatituting t = ee ye 9-2) = malay ae and 


Simplifying yields: 


P| - ¥-+V(1/n+ 1/n)s,° = My- Uys X- Y+ t VU/m+ i/adsy |-¥ 
mei l,u) = (x-F-+tVU/n+ t/n)s,* , x-yr t (1/a+1/n)s,°), 


eet, = x) * Sy. -¥)° 


where Sy (tc 5) 


mee C.1, tor Normal Lly- Ky with oF and or Known 
Peete Z2: pay. bZ25] 


Assuming X and Y are independent and using 


mee Y) - (ily -Ly) 


5 which is distributed N(0,1), as 
V Ox /n + Oy /m 


? 


mie DPivotal-quantity, a 100 7% Gy aI SLCOnStructed. ehus- 


nS 





m= Y) - (i, -jt,) 


Voy2/n + Oy2/n 


Pl 44+Y)/2 < = 7(1+Y)/2! =7. 


Pupstituting z = 21+ )/2 = sta aleay aan and simplifying we 


have: 
mee 3 re 2 =~ =,,.f-2 Daal. 
eR OC, Jn Oy /m = Uy- Uys X-Y+2 Oy /n+ Oy jaler 
SS a ic - _ 2 2 
Sreti,u) = (x-y-2VO,"/nt+ Gy/m, x- yt 2V Oy [aOR a ae 


pyeeec.1. tor Normal C7 / o [Ref. 4: p. 464] 


Assuming X and Y are independent and using 
iB? / me) /(S,"/ a). which is distributed F(m-1,n-1), as 


the pivotal-quantity, a 100 /3 Came Ss eCOMNStrUcted thus: 


2 2 
S/T y 


< F (m-1,n-1)|=‘/. 
Sg. (1+ 7 )/2 


Pp Foy y)/eim-tini) = 


; . i 
pubstituting ewe? = Fasyj/ge bi , and 


Simplifying we have: 


2 
S 8 Z Ox 


x < -——— < va z : : : 
Fas yam tmd) 7 (Sy /Sy JF eiey/2'™ 1,n-1)/=7V. 
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oe 
Sy" /Sy 


Z 2 
or Tew) = (Soy yyaktatemi) (sy [Sy Ps )/gim-tn-t) ). 


eee... £0r Exponential Ax/ Ay = Ly/ Ly Ren an 
p. 466] 
Assuming X and Y are independent and using X A,/Y dy» 
Which is distributed F(2n,2m), as the pivotal-quantity, a 


100 Y% C.I. is constructed thus: 


PI ey y/a6 22m) <— K dy/¥ Ay < Psy y/2(2n2m) | Ses 
Simplifying we have: 


iP Kz/RF +7 ) /262n, 2m) Se <= (Y/RF Oey y/o l2n an) |= Y. 


or (Lrud= (CY/R)F(, yy j2(2n,2m), (YAOE ary y/o(2ns2m)). 


ee tacORY FOR ONE-POPULATION HYPOTHESIS TESTS 

Three types of hypothesis tests are performed in Program 3; 
these are upper-tailed, lower-tailed and two-tailed tests. The 
Beneral procedure used for all three tests 1S the same. 
Basically, a test statistic T, which has an assumed distribu- 
tion, is computed with user-supplied data and then compared 
Oa critical region defined by TS: lass we value is deter- 
mined by the type of test, the assumed distribution, and the 
user-supplied ( value. For two-tailed tests there are 


actually two values of ie used; these values will be denoted 


ae 





here by Ty and T,. Meeomeonrmineg the test then, uf the test 
Statistic T is outside the critical region, then we accept 
Sareoriginal hypothesis H3 otherwise, we reject ee and 
@ccept the alternative H,. 

jemtollowing graphs illustrate the three types of tests 
discussed. These graphs represent a probability density 
function, f(T), where the shaded area under each curve 


represents values from a corresponding cumulative distribution 


femetion, F(T). 


Upper-Tailed Test 


mT) 





il 


it 
O 
A Accept Ao —_——+|.— Reject i: ime 


(Critical Region) 
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Mmemer—lailed Test 
aT ) 





T T 
O 
Re ject H+}. Accept. HH ——_____—____»> 


Two-Tailed Test for Symmetric Distributions (Normal, +t ) 


mcr) 






F(T, )=Q /2 pet) 1/2 


T, il 


i 
i. 
Reject Hope Accept H) —}+ Reject Le 







Median 


Een = 72 
1 - F(T, ) SOL 02 


a al i 


il 2 
t- Reject HT Accept H, ~wtt Reject | 
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Pamupper-tailed test is performed whenever the user 
enters a +1, choosing the upper-tail alternate hypothesis, 
during step one of the solution. A lower-tailed test requires 
tne user to enter a -1 during the solution, and a two-tailed 
Meee requires a 0 entry. For two-tailed unsymmetric tests 
the relationship between the test statistic and the median 
determines which tail 1s actually uSed for the test. If 
the test statistic 1S greater than the median then the two- 
tailed test is performed using the upper tail tr. Value. For 
test statistic values less than the median the test uses the 
lower tail ie value. However, there is no general agreement 
that a two-tailed test should be performed this way. 

The test statistics used in Program 3 will now be listed 
in the same order in which they appear in the User Guide. 
This listing will also include a reference where more informa- 
meenecOncerning a particular subject can be found. Additionally, 
the assumed distribution for each test statistic will appear 
meen that test statistic. 

1. Test Statistic for Normal JL, with Q* Known 


(Ref. 1: p. 431] 


T , using N(0,1) 


(x - Uo dyn 
ee ae 


Meee lest Statistic for Normal Le with a Unknown 


Peet el ep. 451 | 


a0 





he = Sian em adel me @ Tides Sea mia 


Cane 
Ss 9 
X 
Pil eo ne — 250 
meee test Statistic for Bernoulli an eis cee tte) teat 


Ror n > 30 


cl Jo) Se) 


aa lg a N(0,1) 
ey. )s, using 


ememe< 50 the Binomial Distribution in the Statistics 


eo 


Module is used to directly calculate the appropriate probability 
Perecomparison with CY. 
meee lest Statistic for Normal pa with UL Known 
Meret: |p. 452; Ret. 2: p. 104] 


Hor n < 65 


ae ee 


O 


her n > 65 


ee Genre ent cOmy 


on 


5S. Test Statistic for Normal og, with ld Unknown 
Peete il: p. 432; Ref. 2, p. 104] 


PO) < oO4 
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a | using X*(n-1) 
For n > 64 


(pee 2 SG. . using N(0,1) 


2(n-1) 


fest Statistic for Exponential he = ee (eto 


reece 79) | 


nom nm < 52 


AE i ; using aor 


hOr Nn > 32 


r. ae using N(0,1) 


fa 


iWecmeematistic for Poisson oe [Rete 2:e0p,. 248] 


BOT ames 350 


T = n), ; using ¥* (2nx) 





ee using N(0,1) 


See HeORY FOR TWO-POPULATION HYPOTHESIS TESTS 

The two-population hypothesis tests are performed in the 
Same manner as the one-population tests. For an explanation 
of these tests, refer to C above. The test Statistics used 
in Program 4 will now be listed in the same order in which 
Bieyeappear in the User Guide. This listing will include 
the assumed distribution and applicable references. 


fmeetest otatistic for Bernoulli P, 


a. Py ef eee Oe 9 | 


Cay) 


oes nwa iN" using N(0,1) 
a eee ey) Gh) 








( 


meee lest Statistic for Normal Ly = LH, POme a eead Led 





Peete 25> pe 121] 


T = d(vn/s 4), Mommie 1)) tOG m < SL. and 
NtCyietor sie 31 


CA 
Ut 





Bee fest Statistic for Normal Ly = by Ot. amd 
d 
Independent with Oy = Oy pkeee ol: pp. 454-435; 
heme, 23 Dp. 116 


x -y (n-1)sy° + (m-1)sy* 
T = , where s_ = \/—————+*_—____————_ 
ie ol. Pp Lela aa 
S + 
De) ae 


Mims UVC t( ntl -2) ror Wn <= 32. and 


Nee le) ron amie SZ 


meee lest Statistic for Normal Uy = py fOr —X -ancdeet 


: 2 Z 
Independent with Oy ea Oy MiewnmmbRer.. 24 a emeo] 


, using N(0,1) 





meeeiest otatistic for Normal ly = [ly HOE Xecamd. Y 


Se ee ee ere 
Independent with O* # Cy Pvemem 2: 7 a 1419 | 
T = 5 5 mUcMicm CC ecor ds = 50, and 
/s Ss 
+* + Timi rOErdi <i o0. where 
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Ee 


df = largest integer in 








Yn m 


Meee lest Statistic for Normal a, = ar Lom amd. ¥ 


imaependene [Ref. 2: ~ p. 111] 
T = —— msde er (n=) m—1) 


fee test Statistic for Normal P = 0 [Ref. 1: pp. 492- 
eee eee 2 ee. 200] 


T te ee usin wet noo etor 1 < 28, and 


CO oor Nn 7.28 


pee fest Statistic for Exponential IS a ee Ore eta 
imdependent [Ref. 4: p. 310] 


wovicmmo F (ZR 2m ) 


rs 
T 
<i |e 


THEORY FOR METHODS OF APPROXIMATION 


The methods used to approximate distribution values for 


this package of programs come from many sources, and any 


complete discussion of the theory involved would be beyond 


the scope of this paper. Therefore, while every method which 





has been used will be referenced, only those methods unique 
to this programming effort will be discussed in theoretical 
Seeatl, With this in mind, a combination 1isting and dis- 
cussion will follow which will trace the approximation methods 
Meeoein this paper. The approximations unique to this effort 
are the inverse cumulative distribution functions (CDF) for 
the Chi-Square, Student's t and F distributions as well as 
fieeMuitinomial approximation. The inverse Normal CDF 
approximation discussed here was used in Professor Zhena's 
meerams., All other approximations used in this paper are 
Meartacations of methods from the TI-59 Statistics Module. 
For a better discussion of the exact methods used in these 
@iees, the T1-59 Applied Statistics Manual Should be consulted 
Beet. 7}. 

There are two types of inverse CDF approximations used in 
Mmemsepaper. The first two programs use a less accurate 
Type I approximation, while the distribution programs use 
Memmoscliy related, but more accurate, Type II approximation. 
What follows is a discussion of the methods used for 
approximating the inverse Normal, Chi-Square, Student's t and 
meaistrubitons, and the Multinomial distribution. 

ioe Iimverse Normal CDF Approximation [Ref. 5: p. 933] 

a. Type I Approximation 
Iteswapproximation as used in Programs land 2 

for inverse Normal CDF values and as a subroutine for the Chi- 
Square, t and F approximations in those same programs. The 


Type I approximation uses the following set of equations and 
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momotantS tO approximate inverse Normal CDF values, given the 
mgmt probability p. Programs 1 and 2 have the limitation 


fiat p be greater than .5. 


Co + c,t ote oe 5 
Mase atts a mineresm t= 9 Vin(1/(i-p)~) . and 
iE Zz 3 
Cg = Zeng Sila) dy = 1.432788 
c, = oO oo d, = Ils S$) 4le 8, 
C5 = a) OSs 2a5 d. =e) ele GS 


b. Type II Approximation 
The Type II approximation is almost identical 
Memtne, type I approximation. The only difference is the 
meron Of a function which removes the limitation that p 
be greater than .5. This function uses the symmetric 
quality of the Normal distribution and returns the negative 
of the approximation for 1-p, whenever p is less than .5S. 
2. Inverse Chi-Square CDF Approximation [Ref. 5: p. 941] 
a. Type I Approximation 
This approximation is used in Program 1 only 
and uses the inverse Normal CDF approximation described 
above. The inverse Chi-Square CDF approximation uses the 
following set of equations and calculates both Xp) @) 
and cp) (¥) given the input probability p. The reference 
listed above limits degrees of freedom, v, for this approxima- 
tlon to values above 30; however, as can be seen in Appendix B 


Values well below v = 30 produce acceptable approximations. 


ou 





ey) = ieee + Z a)> 
Pp P 


where a 


? 


LO to 
<| 


x 


bee type If Approximation 

liceomappuoOximation 1s used Only in Program 7 and 
incorporates the Type I inverse Chi-Square CDF approximation 
Mmitiean Accuracy Enhancing Technique (AET) which requires the 
highly accurate forward Chi-Square CDF approximation contained 
maetiat Same program. This AET involves taking the output of 
meemtype | approximation and using it in the forward CDF 
approximation to obtain an eStimate, 6, of the actual proba- 
bility achieved by the Type I approximation. This estimate 
1s then used to correct the inverse approximation input, p, 
for any difference between desired and actual probability. 
The corrected inverse input, p', iS computed using the follow- 


me rormula: 


pep sep =p) 


femuoing p' as the new input for the inverse CDF approxima- 
tion, a more accurate approximation is achieved (see Appendix B). 
This AET is not used for v > 30 since the Type I approximation 
Meethen quite accurate. 
See Inverse Student's t CDF Approximation 
a. Type I Approximation 
Mmcomanpneximation 1S used in Programs i, 2, 8 
and 9 and requires the inverse Normal CDF approximation 


already described. The following equation represents one of 
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eral approximations to the inverse t developed by Professor 
Donald P. Gaver, Naval Postgraduate School. This approxima- 


Sees limited to values of v > 2 [Ref. 8]. 


b. Type II Approximation 
iis eapeGoxiMatilon 1S used Only in Program 8 
and, like the Type II Chi-Square approximation, uses an AET 
tO increase the accuracy of the Type I approximation. The 
Meeeedires for this AET are exactly the same as in the Chi- 
Be@eare approximation. Additionally, the restriction on v 
has been removed by using the following relationship to 
generate inverse approximations when v = l. 
t,@) = PSRTISST , Where for p < .5, 2p is used in place of p 
fee inverse F CDF Approximation [Ref. 5: p. 947] 
@. tfype i Approximation 
(icmp reOximatvon 2s used in Programs 1 and 2 
and uses the inverse Normal approximation discussed earlier. 
The following set of equations is used to generate the 


approximation: 





fee 2 
2w ee  * 


FL (vy + Vp) =e weviere Ww = , ke 


i 


bee lype il Approximation 
Uiisedpproximation uses the ABT discussed earlier 
along with two more techniques to provide inverse approxima- 


meets £LOr all values of V4 and V5° For the case where v., or 


1 
a = 1, the Type I inverse t approximation discussed earlier 


aS used in the following relationships to generate the required 


inverse F CDF approximations: 


FL (,v) end ies , or equivalently 


(1+p)/2 


me(v,i) = ee 
Dp ’ teepy/o'” 


For the case where Nts 1, the following relationship 


is used: 


it 
Fo(tst) = tantimp) - 
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S. Multinomial Approximation 
The method used to generate multinomial density values 
is not truly unique; however, it uses the following equation 
in a way that minimizes rounding errors: 


N! 
ey; 5, ---sn,) = Th, PV (ngh)- tn, t) (Py 2) (Py°2)., accra) 


Computations are accomplished in the following order to 
avoid, as much as possible, multiplying extremely small 


Values by extremely large values: 


mee 1 )\ /(N - 1)(p."2) N-k Dy 
1 shila ala cae neat ad Cleese! 
a 2k Dy 
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ive 11-59 





iifeschoice of the TI1-59 as the calculator for this pro- 
cramming effort was based on two factors. Professor Zehna's 
Original programs were written for the TI-59; and each student 
in the Qperations Research (OR) curriculum is tssued a TI-59 
for use in basic probability and statistics courses. In 
general, the use of hand-held programmable calculators has 
been shown [Ref. 6: p. 1] to increase student learning and 
Capability. Further, it is intended that the programs 
described in this paper will be used by OR students in their 
coursework. 

Using the TI-59 offered some disadvantages and some 
advantages in developing the programming package presented 
here. The disadvantages of limited storage, slow computation 
time and awkward data entry sequence are discussed elsewhere 
in this paper. There are two advantages of the TI-59, how- 
ever, that deserve noting here. The programming steps and 
procedures used in the TI-59 are easily learned and logical. 
This ease of programming makes complex computational methods 
easy to program. Another advantage of the TI-59 is its 
compatibility with the PC-100C printer. Using the printer/ 
Bereulator combination greatly simplifies writing, editing, 
and error diagnosis. These advantages of the TI-59 make 
programming relatively easy and also allow the capabilities 


of the TI-59 to be used more fully. 
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V. ALTERNATE SOLUTIONS 


ies problems presented in this paper could have been 
solved in any number of equally valid ways. This section 
Smemeprietly discuss the specific alternatives which could 
Meerease the capability of the solutions employed here. 
This discussion will first focus on alternative hardware 
Which could be used, and then on program changes which might 


be made. 


A. HARDWARE ALTERNATIVES 

ms discussed earlier, the use of the TI-59 calculator is 
appropriate for this programming package; however, the methods 
and techniques used in this paper are equally suited to other 
programmable calculators or computers. Indeed, the use of 
any Other calculator or computer with more storage capability 
than the TI-59 might be a better vehicle for the package of 
programs presented here. These programs require a total of 
14 magnetic cards (22 separate sides) using the TI-59. 
Current console model computers could easily store all of 
these programs at the same time. 

The TI-59 has the ability to use program modules, and 
the nine programs from this package could easily be combined 
to form the basis of a new module. By careful elimination of 
redundant functions, such a module could also accommodate an 
ANOVA and/or regression package. A module of this type might 


find widespread acceptance, especially in the classroom. 
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B. PROGRAMMING ALTERNATIVES 

The choice of distribution approximating methods used 
here and the general layout of the entire package might 
both be improved with some additional effort. The methods 
of approximation used here were chosen rather arbitrarily. 
it is possible that more appropriate methods of approximation 
exist. More appropriate methods might include methods with 
fewer program steps and equal accuracy, as well as methods 
Which take less computation time with equal accuracy. 

The current layout of the nine programs in this package 
might be improved by reducing the number of programs. The 
present program size precludes this; however, by eliminating 
Some program functions, it might be possible to organize all 
of the inverse CDF approximations in one program. Shorter 
approximating methods might make possible a similar program, 
Sutewith both inverse CDF and regular CDF approximations. 

The extensive user guides with their sample problems are 
useful for students, but could prove awkward for a more 
experienced user. An obvious alternative would be a shortened 


user guide directed at those more familiar with the programs. 
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Vie CONCEUSION 


Two problems were presented for solution in this paper. 
Qne problem involved expanding and modifying an existing set 
of TI-59 programs into a user-friendly package. A second 
problem involved developing a set of distribution approximating 
programs. The solution to the first problem incorporated 
increased capability, standardized data entry, and detailed 
user guides into a package of nine programs. The first 
priority in this solution was providing a format compatible 
With a student's needs while maintaining the capability 
mBequired by a more experienced user. The TI-59 calculator 
proved to be a very useful tool in this solution, and 
demonstrated the generally unused capability of the current 
generation of hand-held programmable calculators. 

The approximation programs presented as a solution to the 
second problem mentioned above provide accurate and compre- 
hensive approximations. These programs practically eliminate 
the need for tables of values and solve the interpolation 
problem present in all such tables. 

Together these two solutions represent a package with 
considerable capability in computing confidence intervals, 
performing hypothesis tests, and generating approximate distri- 
bution values. A comprehensive set of user guides makes 
this same capability available even to inexperienced users. 


The methods used in preparing this TI-59 programming 
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eee Cerlyecappltcable to other calculators or 
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APPENDIX A 
PROGRAM 1 USER GUIDE - One-Population Confidence Intervals 
INTRODUCTION: The purpose of this program is to compute 100°Y% 
confidence intervals (1,u) or bounds fa and y| for the following one- 


population situations: 


NORMAL U with CO known 
NORMAL [LL with O* unknown 
BERNOULLI op 

NORMAL O* with LL known 
NORMAL O* with [U unknown 


EXPONENTIAL vor LL . 

The routines in this program require percentiles from either the 
Normal, Chi-Square, Student's t or F distributions. Bach routine will 
automatically generate an approximate percentile; however, when 
additional accuracy is desired or small sample sizes are involved the 
use of percentile values from either standard tables or the distribution 
approximating programs is recommended. In step two of each routine the 
user can choose to accept the approximate percentile, by storing the 
appropriate percentage in storage register 09 (Rog) or he can store 
the percentile value in Rag: Some routines also require percentile 


values in Ry 3° 
GENERAL PROCEDURES: 

1. Use any library module, and after reading all three card sides, 
press [ D | comme Datta ron s(o39 .39)) . 

2. For data entry press [ p | followed by data point Xs » Xs 


R/S], etc. for each x; (i=1,2,...,n) until all points have been 


WT 





entered. Mistakes in data entry should be corrected immediately by reen- 


tering the unwanted point and pressing E+], then enter the 
correct data point and press . Alternate data entry using 


summary statistics is detailed in applicable routines. 

3. For one-sided confidence bounds rather than intervals replace 
(1+°7)/2 with Y and (1-))/2 with 1-Y everywhere they appear 
e.g» 2147 )/2 becomes aye) and proceed as usual, ignoring 1 or u 
as appropriate. 

4, When solving consecutive problems, care should be taken to clear 
all previously used registers. Pressing | D | Will clear all registers. 
PROGRAM 1 SPECIFIC PROCEDURES: 

G.I. For NORMAL [J with (/* known 


1. Enter data using eee aC esiuOne | smerny Alternate entry: 


07 | 


store n in Ro3) x in Rog anges ik 


0 07) 
2. Store either (1+))/2 in Rog OR (14 )/2 in Ry, 


3. Press 1 is displayed, then 
press u is displayed 


G.I. For NORMAL [{ with O* unknown 


1. Enter data using [ D |(Alternate entry; store n in Ro3: 


s inR and x in Rog) 


07 
2. Store either (1+7)/2 in Rog OR bey feet) in Ry, 


3. Press 1 is displayed, then 
press u is displayed 
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See ror BERNOULLI p 


1 a 


a. 


. 


Seore 9nx in Roy and mn 17 Ro3 


Store either (14+7)/2 in Roo OR 9 (Ge? 2 = aaa a in Ry, 


AND 


Wgeavgy/2\ctee sce) in Ry 3 
Press 1 is displayed, then 
press uis displayed 


G.I. For NORMAL G* with (f known 


a 


Enter data using [ D | and store iO ra Rog (Alternate entry: 
ear ; 
store Dace in Ryo, n in Ro3 and LL in Rog) 
~~ : : Z : 
Store either (1+/)/2 in Rog OR (a-yy/2 in Ry, 
AND 
Z (n) in R 
X (1+y)/2 13 
Press 1 is displayed, then 


press u is displayed 


G.I. For NORMAL * with unknown 


1. 


Enter data using | D | (Alternate entry: store n in R93 and 


fet )s“ in Ryo) 
Store either (1+Y)/2 in Rog OR X-yyyeleh ine 


AND 


2 _ . 
Press 1 is displayed, then 
press u is displayed 
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%& 
C.1. For EXPONENTIAL jor LL 
1. Enter data using| D | (Alternate entry: store n in Ro3 and 
rarer Rog) 
2. Store either (1+ )Y)/2 in R,g OR (2n) in R 
SeLNeH 09 SX a-y)/2 11 


AND 


X ey /2(2") in Ry3 
3. Press 1 is displayed, then 
press u is displayed 
* 4, To compute a confidence interval or bound for eles 1/)\ press 
[ B'| rather than above. 
PROGRAM 1 ADDITIONAL CAPABILITIES: 


Inverse Normal CDF Approximation 





im otore p in Rog (ee a) 


2. Press a is displayed, and 


Z is available in R 
(Sp) 


13 


Inverse Chi-Square CDF Approximation 





i Store p in Rog (pageant \yeand store vin Ray Gest) 


2. Press en is displayed, and 


2 ; ; 
De Gea is available in R,, 


Inverse F CDF Approximation 





meeestore p in Rog (p > .5), v, in R,, and eee (v,71,vof1) 


2. Press [RCL] F (vy 1Vp) is displayed 


Inverse Student's t CDF Approximation 





ie Store p in Rog (ei eo aand vin Rig ca 


2. Press tiv) is displayed 
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PROGRAM 1 


A 


B 


PROGRAM 
OO 
O1 
O02 
03 
O4 
05 
06 
07 
08 
09 
10 
Leh 
ie 


1 
14 


LABELS USED: 


1 STORAGE REGISTER CONTENTS: 


clear 


Ee 
n 
used 


clear 

evear 

Emon GC). 

x or LL 
ee) )/ 2) or p 
used 

CDF value 
aise 

CDF value 


Vv 


LS 
16 
me 
18 
19 
20 
Al 
22 
23 
24 
25 
26 
27 
28 


29-49 


Si 


RCL 


SLO 


used 


clear 


clear 


clear 


clear 


clear 


used 


clear 


clear 


used 


used 


clear 





PROGRAM 1 SAMPLE PROBLEMS: 
1. Suppose x= 69.7, n= 8, O* = 3.5. Find a 90% C.1. for 
estimating the mean. 
SOLUTION: 
a) Store n=8 in Ros: = 69e7 in Rog and OG = Broan 
(2) Store (1+/)/2 = .95 in Rog OR 2 ge = 1.645 in R,, 


ee ftess then to display 


1 = 68.61179492 1 = 68.61193477 


07 


OR 
70 .78820508 as u = 70.78806 523 


u 
2. Given the following observations from a population with a known 
variance of 15, find a lower 99% confidence bound on the mean. 
mopet7S 160 199 167 145 157 182 192 165 
SOLUTION: 
(1) Enter data using [ D |sequence and store (= “fi5 in Ron 
(2) Store “Y= .99 in Rog OR Zigg = 2-326 in R 


(3) Press then to display 


i=l Oomieezolo i, Ss alieys! ot Sry tig 
OR 
ignore u a ignore u 


iL 


3. Find a 90% C.1I. for the mean of the following test scores. 
35 34 46 20 38 38339 32: 49 41 25 «1843 51 38 42 «6029 59 53 27 33 
SOLUTION: 
(1) Enter data using | D | sequence 
(2) Store (1+ ))/2 = .95 in Rog OR t go(19) = 1.729 in R,, 
(3) Press then to display 
1 = 3339962583 eee 7 ae 
41 80037417 a u = 41.80609048 


i 


Ik 


u 
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4, Find a 95% upper confidence bound for fit given the following 


information: n= 55, s=i5 and x = 85. 


SOLUTION: 
OD) Store n = 55 in Ros: s=15in Roo and x = 65 in Rog 
(2) Store “Y= .95 in Rog OR t 95(#) = 1.673 in R,, 
(3) Press then to display 
ignore 1 ignore 1 
OR 
u = 88.38398464 l=see-cejo0on! 


5. Of 1000 people treated with a certain drug 200 showed a reaction. 


Find a 90% C.1I. for the proportion of the sample population that 


Will show a reaction. 


SOLUTION: 
(1) Press | D | then store nx = 200 in Ry, and n = 1000 in Ry, 
(2) Store (1+Y)/2 = .95 in Rog OR F 95(402,1600) =1.13 in R,, 
AND 
F 96(1602,400) = 1.14 in Ry, 
(3) Press then to display 


1 


.1793622306 i 
OR 
-2219 574577 u 


.1796719191 
fe O25 


Uu 


oS 





eo. 


Find a 95% upper confidence bound for the proportion of the sample 


population if nine of 24 treated were affected. 


SOLUTION: 
(1) Press oa then store nx = 9 in Ro4 and n = 24 in Ro3 
(2) Store Y= .95 in Rog OR F g5(20,30) = 1.93 in R,, 
AND 
F 95(32,18) ete eat) 343 
(3) Press then to display 
ignore l ignore 1 
OR 
u = .5629409878 u = . 5626822157 


feenc, a 95% C.1. for Oka and for (Q given the following observations 
from a sample population with [L= 65.0 . 
100 15 73 46 65 98 79 38 68 85 
SOLUTION: 
(1) Enter data using [> | and store [{= 65 in Rog 
= : Z 
(2) Store (147 )/2 = .975 in Rog OR X 1025610) 
AND 
2 1 = i 
x 1975610) 20.48 in Ry 3 


(3) Press then to display (for o*) 


428 .7608301 1 = 428.8574219 
OR 
2726 .1 24604 . u 


3525-10 


aa 


Ea! 
it 


& 
tt 


2702 461 538 


Taking square roots the following limits for GO are displayed 


J 
I 


20.706 54076 1 
OR 
52 .21230319 u 


20 . 70887302 


or 
| 


bimge52d> 


Sy 





S. 


9. 


Suppose n= 15, prac = 88476 and LU = 30.5. Find a 97.5% upper 
confidence bound on the standard deviation. 
SOLUTION: 
(1) Press [ > | then store peas =" Gc476 in Roo» n= 15 in R93 
and fe Bono in Rog 
= : 2 
(2) Store “Y= .975 in Rog OR xX 1025619) 
AND 
2 ie : 
iG 975615) 27 49 in Ry3 


(3) Press Tien and take square root to display 


ignore 1 ignore l 


6.26 in R 


i 


OR 
u = 109 .2668647 —— u = 109.1078035 


Find a 95% C.I. for JO given the following observations. 
mo 15 73 46 65 98 79 38 68 85 
SOLUTION: 
(1) Enter data using oe 
(2) Store (1+Y)/2 = .975 in Rog OR X* 990(9) eae 
AND 
2 _ ; 
4 97569) = 19.02 in Ry, 


(3) Press then and take square roots to display 


1 = 18.3742 1 = 18.54896609 
OR 
4947192697 u 


I 
II 


I 
I 


HO G231 501 5k 


u 
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10. Suppose n = 45 and s* = 36 find a 90% lower confidence bound for 


Oe. 
SOLUTION: 
(1) Press then store n = 45 in Ro3 and (n-1)s~ = 1584 
in Rao 
(2) Store ‘°Y = 90 in Rog OR X* 45 () = 32.5 in R,, 
AND 
2 ' > 1 2 
x 90) = 56.4 in Ry 3 
(3) Press to display 
1 = 28.10395538 1 = 28.08510638 
O 
ignore u =i Lenore 


11. Given the time tc failure of an electron tube is an exponential 


random variable, and the sum of 25 times to failure is 25242. 


Find a 95% C.1. for ). 


SOLUTION : 


(1) 


(2) 


(3) 


Press El then store n = 25 in R,. and x = 25242/25 


03 
= 4 2 = 3] j 
Store (14° 7 )/2 = .975 in Rog OR X 10256) = 31.92 in Roy 


AND 


X* .975(0) = 70.92 in Ry, 


Press then to display 


1 = .0006407042 1 = .0006322795 
OR 
.001414878 _ u 


u 001404801 5 


26. 





12. Six expensive pieces of equipment had the following times to 
failure. Find a 95% C.I. for the mean time to failure. 
Peero 3119.0 258.3 1402.7 612.9 2211.2 
SOLUTION : 


(1) Enter data using [> | 
(2) Store (1+ ))/2 = .975 in Rog OR Y* y96(t2) 
AND 


XC opel) = 23.34 in R,5 


(3) Press then to display 


ees O71 6219738 i =" 671 6109683 
OR 
3578 .128804 — u 


4.40 in R 


el 


3.562 . 590909 


u 


oa 





PROGRAM 2 USER GUIDE - Two-Population Confidence intervals: 
INTRODUCTION: The purpose of this program is to compute 100°Y% 


confidence intervals (1,u) or bounds E and u| for the following two- 
population situations. 


BERNOULLI Py-Py for large sample sizes 


= : : 2 
NORMAL [Ly [ly for X and Y paired with Q~ unknown 

: : ae ape ee 2 
NORMAL [ly Hy for X and Y independent with G,"=Gy ei 

2 : : 2 Oe 
NORMAL Ly py for X and Y indep with Ox and Cy Known 
NORMAL Ore Oy" Eon Aen dumectncependecnit 


EXPONENTIAL dy/ Ay = fy/[L,y for X and ¥ independent 

The routines in this program require percentiles from either the 
Normal, Student's tor F distributions. Each routine will automat- 
ically generate an approximate percentile; however, when additional 
accuracy is desired or small sample sizes are involved the use of 
percentile values from either standard tables or the distribution 
approximating programs is recommended. In step two of each routine the 
user can choose to accept the approximate percentile, by storing the 
appropriate percentage in Rog: or he can store the the percentile value 
in Raa: Some routines also require percentile values in R43" 
GENERAL PROCEDURES : 

1. Use any library module, and after reading all three card sides, 
press [p | to repartition (719.29). 

2. Data entry for independent data (DEI Sequence), press [ > | 
followed by data point Xs [z/sl, Xe ays , etc. for each x, (ads 2; 


+ee,n). When all of the x, 'S have been entered press followed by 


data point Yas [R/sl., Yeua? [z/s] , etc. for each 5 (i=1,2,...,m). 
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When all of the y,'s have been entered press [cro]. Mistakes in 


data entry should be corrected immediately by reentering the unwanted 


data point and pressing [E+], then enter the correct data 
point and press . Alternate data entry using summary 


statistics is detailed in applicable routines. 

3. Data entry for paired data (DEP Sequence), press| D |followed 
by a? [xxt] , Ys , eee [at], Ss4q? R/S], etc. ior each data 
pair Xe Vs (4=1,,2, ; ees, then press » Mistakes and 


alternate entry are as above. 

4. For confidence bounds rather than intervals, replace (147 )/2 
with “Y everywhere it appears (e.g. tay 2”? becomes a ) and 
proceed as usual ignoring either 1 or u as appropriate. 

5. When solving consecutive problems, care should be taken to 
clear all previously used storage registers. Pressing | D | Will clear 


all registers. 


PROGRAM 2 SPECIFIC PROCEDURES: 


C.1. For BERNOULLI Py-Py for large sample sizes” 


1. Press ioe then store nx in Roa? my in Roo ea Roa) m in Ron 


2. Store either (1+))/2 in Rog OR 21) /2 in Ry4 


pee Press 1 is displayed, then 


press |xzt}| u is displayed 


* Large here means NPys MPys n(1-p,), n(1~-p,) all greater than five. 


oo 





Cee ror NORMAL [Ly Ly for paired X and ¥ (n pairs) 


ile 


4, 


inter data using DEP Sequence, degrees of freedom, v = n-1, will 
be displayed. (Alternate entry: store X-y in Ryo? s/s) 2 in Ryo 
and n-i in Ry 9) 

Store either (1+))/2 in Rog OB bey ysete-2) in Ry, 
Press 1 is displayed, then 


press |xzt{/ u is displayed 


When differences, i are given use Program i. 


G.I. For NORMAL {U,-[l, with g,°= gy” = @* unknowm 


is 


Enter data using DEI Sequence, degrees of freedom, v = n+m-2 is 


displayed. (Alternate Cllviay GfamencOome iM itl Ree, Mm in Rin X-y sia) 


t> 03 
R ye a) in R y(y -y)* in R,, and (n+m-2) in R,,) 
Or zi 16 i 26 “19 

: : on 
(Alternate entry #2: Store n in Ro) ee in Roy? ey in Roo: 

: : eis 

then press| D'], store m in Roa? we ei Rox and oe 2 Roo: 
then press , v is displayed) 
Store either (1+)/2 in Rog OB tay) /otata2) in Ry, 


Press 1 is displayed, then 
press u is displayed 


C.l. For NORMAL [ly- Ly with ore and aa known 


1. 


Enter data using DEI Sequence, tnen store anys oer) Ran and 


C,F/m in Ron (Alternate entry; store x-y in Ryo and 


V\o,2/a) + (G,2/m) in Rp) 


puene Clumer Gia) /e atid Rog OR aides) /2 in Ry4 
Press 1 is displayed, then 
press u is displayed 


6.0 





G.I. For NORMAL g,°/ gy” 


1. Enter data using DEI Sequence (Alternate entry: store n in R 


a. ; ae 
sy, in R,o, m in Ro, and s,~ in Roo) 


2. Store either (1+7)/2 in Rog OR Bray 2a) in Ry, 


oy 


AND 


aye) / 2h yn-1) et 13 
3. Press 1 is displayed, then 


press |xzt/ u is displayed 
G.I. For EXPONENTIAL A,/)y = Lly/[Ly 
1. Enter data using DEI Sequence (Alternate entry: store n in Ry os 
x in R,g, m in Ro3 cand Rog) 
2. Store either ea) (e in Rog OR SCrayay ena) in Ry, 
AND 
a aye) in Ry 4 


3. Press lL is displayed, then 


press |x2t/ u is displayed 
PROGRAM 2 ADDITIONAL CAPABILITIES: 


Inverse Normal CDF Approximation 





Mm store p in Rog (p => .5) 


@. Press to display 7 


Inverse Student's t CDF Approximation 





1. Store p in Rog (pet end vy in Rig Cae) 


2. Press to display ey) 


Inverse F CDF Approximation 





ime store p in Rog (eee 25); v, in Roy and v, in Ry, (v,#1,V#1) 


2. Press to display Fp(v1,v2) 


B. 





PROGRAM 2 LABELS USED 


A A’ GTO BRE 
B Se RST 
C = SBR 
D D' STO 
B E RCL 


PROGRAM 2 STORAGE REGISTER CONTENTS: 


OO clear ion 

O1 nx or Ly; 16 used 

02 my or enc 17 aS or Gare 
03 norm See 

O4 mor eg 19 vor (n-1) 
i) x,“ | 20 v, 

em) x,y, Zia 

07 ee or C,£/m 22 clear 

o8 y 23 used 
Beep or (1+ 7)/2 24 used 

10 x-y 25 clear 

11 CDF value 26 clear 

2 s lf 27 used 

13. CDF value 28 used 

14 used 29 used 
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PROGRAM 2 SAMPLE PROBLEMS: 


1. 


In a survey of 400 people from one city 188 preferred Brand A soap 
to all others; and in a sample of 500 people from another city 210 
preferred the same product. Find a 95% confidence interval for 
Py Py : 

SOLUTION : 

(1) Press [ D | then store nx = 188 in Roy» my = 210 in Roo, 


n = 400 in Ro3 and m= 500 in Roy 


(2) Store (1+7)/2 = .975 in Rog OR 2 996 = 1.960 in R,, 
(3) Press then omens lan 
1 = -.0153123448 1 = -.0152991877 
OR 
u= .1153123448 = ei COSTE 77 


In a survey of dieting effects, ten women were selected, weighed 
and placed on a diet for two Weeks. At the end of that time they 
were reweighed. The results are listed below. Find a 99% confidence 
interval for the difference in means. Also find a 95% upper 
confidence bound for the difference in means. 
Bepore 119 i122 136 130 129 136 134 133 119 115 
After 114 119 134 126 119 137 124 4127 #119 107 

SOLUTION Part a. 

(1) Enter data using DEP Sequence, 9 is displayed 

(2) Store (1+) )/2 = .995 in Rog OR b 9g5(9) = 3.250 in Ry, 


(3) Press then to display 


1 = 7054594846 1 = .7328694249 
OR 
8.694 540515 —_ u 


1 
I 


Uu 


8 667130575 


en: 





SOLUTION Part b. 
(1) Enter data using DEP Sequence, 9 is displayed 
(2) Store “Y= .95 in Rog OR tg «(9) eee aa 


(3) Press then to display 


ignore l ignore 1 
OR 
u = 6.92957002 u = 6.937461644 
3. Suppose x = 10, y = 5, n= 25, ana S53 = 10; find a 90% C.1. for 
Ly-Uy - 
SOLUTION; 


(1) Press [ D | then store x-y = 5 in Ryo; Sq/a/n = 2 in R,, 
= - 2 j 
and (n-1) = 24 in Rig 


M2 store (1+) )/2 in Rog OR t gn(24) = 1.711 in R 
(3) Press then to display 
ee 50) 603572 eee 76 


OR 
8 418196428 u = 8.422 


ale 


I 


u 


4. Two small classes used different, but equivalent methods on a 
common exam. Their scores are listed below. Find a 90% confidence 
interval and a 90% lower bound for the difference in means. 
es 82 87 91 HH 97 76 & 98 92 57 80 53 &} 

MH 91 6294 92 87) =—979 86 7590 73 «28393 6589 68) =e 
SOLUTION: Part a. 
(1) Enter data using DEI Sequence, v = 27 is displayed 
(2) store (1+7)/2 = .95 in Rog OB & 9 6(27) = 1.703 in R,, 
(3) Press then to display 


l = -12.62699316 Ll = -12.63357801 


5828916241 "= u 


u 


5 835501092 
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SOLUTION Paxt b. 

(1) Enter data using DEI Sequence, v = 27 is displayed 

(2) Store “Y= .90 in Rog 

(3) Press to display 
1 = -10.49428744 1 = -10.52422048 


OR 
Lgnore u ignore u 


OR t 99(27) = 1.314 in Ry, 


Suppose n = 16, m= 15, X= 14.3, y = 12.5, )(x,-x)* = 67.2 and 
(y,-¥)* = 103.7 . Find a 99% upper confidence bound for the 
difference in means. 
SOLUTION: 
cl) Press [ D | then store n = 16 in Ryo m= 15 in Ro3) 
1.8 in Ryo» ¥(x, -x)* O62. Ry 6) Y(y,-¥)* = 103.7 
in Rog and (nt+m-2) = 29 in Rig 
(2) Store “Y= .99 in Rog OR © 9929) = 2.462 in R 


(3) Press [¢] then to display 


ignore 1 ignore l 
OR 
= 3.959764742 u = 3.94800 5095 


x-y = 


ial 


Emppose n = 32, m= 10, x = 13, y = 15, Ges = oan oh = 129} 
find a 95% confidence interval for jly-[Ly: 
SOLUTION: 


fee) Press | D | then store x-y = 2 in Ryo and 


V(0,2/n) + (o,2/n) = f3 in B,, 


(2) Store (1+7)/2= .975inR,, OR 3% = 1.960 inR 


09 975 fee 
(3) Press [ a | then [xt } to display 
SOS SESS) 1 = -1.394819583 
OR 
ee elo a u= 5.394819583 


6.5 


Given the following sets of data with known variance, estimate the 
difference in means with a 90% confidence interval. 
@ = 3-7 and gy = 2.5 
mero 447 «62.04 1.7 6.2 7.2 3.6 
Mme 7.0 4.3 6.5 8.7 4.5 6.1 8.6 
SOLUTION : 
(1) Enter data using DEI Sequence and store Ge oe 7am 
7 and GF /n = 257/18) aba 
(2) Store (1+7)/2 = .95 in Rog 
(3) Press [a] then [xxt] to display 
een (ee Dt iz -3.77112862 
- .7 536774686 - u 


O70 


OR 2 =elOt5 2m) kh 


9D de 


i 
tt 


iD 


- «7538713802 

Find a 99% confidence interval for ¢y,*/Q,* using the following 

data. 

Kk: 82 8791 5497 76 64 98 92 57 80 53 & 

mee ce 94 92 87 8979 86 7590 73 «48393 6589 68 
SOLUTION: 


(1) ter data using DEI Sequence 


(2) Store (1+7)/2 = .995 in Rog OR F 99 5(12,15)=4.25 in Ry. 
AND 
F g95(15,12) = 17240 R43 
(3) Press then to display 
1 = .38646389829 He ewe oe St 0S 
OR 
= 7.818933411 u = 7.773767765 
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9. Find a 95% confidence interval for Gf,"/ gj," given n= 8, m= il, 


Zz Zee 
Sy = 4 and Sy BO 


SOLUTION: 


(1) Press [ >] then store n = 8 in Rao m= ii in Ro3) 


aa oo : 
as 4 in Ryo and ao Sorin Ron 


(2) store (1+Y)/2 = .975 in Rog OR F o9.(7,10) = 3.95 in Ry 


AND 
F g75 (10,7) = 7) an Ry 3 
(3) Press then to display 
1 = .280241493 1 = .2812939522 
OR 
u = 5.35837828 a u = 5.288888889 


10. Find a 90% confidence interval for fue ele given the following 
data; 
fey )6068old«Cc12 «ik 312 «215 
fey 15 4 66 lh 13 
SOLUTION: 


(1) Enter data using DEI Sequence 


(2) Store (1+77)/2 = .95 in Rog OB F gs(12,14) 2.55 in R, 


i 
AND 


F gs(14,12) = 2.64 in R 


(3) Press then to display 


1 = .3688954123 1 = .3822288409 
OR 
u = 2.472430016 | u = 2.573164557 


13 


Ge 





tac. 


Two types of electric bulbs are observed as to length of life, 


yeilding the following results. Are the means significantly 


Sarerent? “(CY = .1) 


Mepet n= 46 x = 1070 
jepe2 m= 62 y= 1041 
SOLUTION: 
(1) Press [>] then store n = 46 in Ryo» m= 64 in Ro3) 
a7 Olen Ryg and y= Ont in Rog 
(2) store (1+7)/2 = .95 in Rog OR F 5.(98,128) = 1.33 inR,, 
AND 
F 95(128,98) =e Oneia Ry 3 
(3) Press then to display 
1 = .7102715868 ley oleOlo sie 
u = 1.345613 48 u = 1.362056075 


Thus With 90% confidence the means are not different. 
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PROGRAM 3 USER GUIDE - One-Population Hypothesis Tests 
INTRODUCTION: The purpose of this program is to perform one-tailed or 


two-tailed hypothesis tests for the following one-population situations: 


NORMAL fee With o* known 
NORMAL Ll, with O* unknown 
BERNOULLI Pp, 

NORMAL on with LL known 
NORMAL CO, with {£ unknown 


EXPONENTIAL jl. = oe 
POISSON eo 
GENERAL PROCEDURES : 

1. This program requires the Applied Statistics Module. After 
reading all three card edges, press to repartition (719.29). 

2. For each new problem press to clear all registers and to 
prepare for the following parameter entries. Enter o. and press [z/s], 
enter CY and press , then enter either 0, -1 or +1 for the desired 
alternate hypothesis, as in the table below, and press [z/s] . 

0 for Hy: 6 # @,, Two-tailed test 


ee tor H,: gO =< Or. Lower-tailed test 


it 
ci, ior Hy GO > Go: Upper-tailed test 
3. For data entry press [ p | followed by data point Xe, ; 
Xd? , eta for each x, (i=1,2,...,n) until all data points have 


been entered. Mistakes in data entry should be corrected immediately by 


reentering the unwanted data point and pressing ener 


enter the correct data point and press . Alternate data entry 


using summary statistics is detailed in each subroutine. 


es. 





4, At the conclusion of each test a 1 for Reject or a OQ for 
Accept is displayed. The significance level for each test is usually in 
the T-register (Rin) § however, for two-tailed unsymmetric tests the test 


Piabrstic stored in R,,. is used as follows: if a O is displayed tnen 


eZ 
the test must have failed to reject in both tails. If a 1 is displayed 


then the test Was upper-tailed if R,. is greater than the median and 


2 


lower-tailed if R is less than the median. 


eZ 
PROGRAM 3 SPECIFIC PROCEDURES: 


Tests for NORMAL jf, with O* known 


1. Press EE enter a press [r/s], enter (, press [R/s] , 
enter 0, -1 or +1 and press ; 


2. Enter data using [> | and store GQ in Rog (Alternate entry: 


store n in Ro3: xin R,. and GO in R 


10 07) 
be) Press either 1 (reject) or O (accept) is displayed, 
press to display significance level. 


Tests for NORMAL pee with Of unknown 
1. Press eae enter eb press [z/s], enter CY, press [z/s] , 
enter 0, -1 or +1 and press ; 


2. Enter data using [ D | (Alternate entry: store n in Bo! eo eT) Rao 
and s, in Ray): 
3. either 1 (reject) or 0 (accept) is displayed, then 


press to display significance level. 


7D 





Tests For Bernoulli Py 


i 


Press Bea, enter Py» press [r/s], enter C{, press [r/s], 
enter 0, -1 or +1 and press ; 


Store nx in R., andninR.. 


Ol 03 
Press either 1 (reject) or O (accept) is displayed, then 


press bemdisplay significance level, 


Tests For NORMAL ne with LL known 


i. 


3. 


Press ere enter Ones press [R/s], envemicy, press [e/s], 
enter 0, -1 or +1 and press ‘ 


Enter data using [D] and store ee cel Ry 9 (Alternate entry: 
, : Eee 
store n in Ross es in Roy er in Roo and fl in Ry) 
Press either 1 (reject) or 0 (accept) is displayed, then 


press to display significance level. 


Tests For NORMAL Ge with LL unknown 


i. 


Press | B'], enter g.* , press [R/s], enter CY, press [R/s]_, 
enter 0, -1 or +1 and press : 


Enter data using[ D | (Alternate entry: store n in Ro3 and 
-1)s * ; 
(n 1)s, in Rio): 


Press either 1 (reject) or O (accept) is displayed, then 


press to display significance level. 


Tests For EXPONENTIAL {l. = tee 


Hee 


Press Para enter {i 9: press [z/s] , enter CX, press [z/s], 
enter 0, -1 or +1 and press ‘ 


Enter data using [ D |(Alternate: store n in Ros and x in Ry) ‘ 


Press either 1 (reject) or 0O (accept) is displayed, then 


press to display significance level. 








Tests For POISSON } . 


1. Press ee enter ees press [R/s], enter CZ, press [R/s], 
enter 0, -i or +i and press : 


Enter data using | D |(Alternate entry: store n in R 


nx in Roy): 


03 and 


Press either 1 (reject) or 0O (accept) is displayed, then 


press to display significance level. 


PROGRAM 3 LABELS USED 


A 


B 


Aa 


Be 


TAN 


SIN 


COS 


CLR 


GRD 


SUM 


CE 


PROGRAM 3 STORAGE REGISTER CONTENTS 


00 
O1 
02 
03 
o4 
05 
06 
07 
08 


09 


Q, -1 or +1 
moe Or eo 
D Z 

xe 

af 

n 
clear 
clear 


degrees of freedom 


CT 


Z N 
or r 
K OG, por A, 


clear 


10 
ah 
eZ 
13 
14 


LNX 


LOG 


Bue 


INT 


poroa [i ox GeD or 
CL 

test statistic 
clear 


= 


15-29 used 


Ee 





PROGRAM 3 SAMPLE PROBLEMS: 


i. 


According to an early encyclopedia the average rainfall in a city is 
30.1 inches. Rainfall during the past five years has been: 
Boo 27-0 35.1 32.60 «25.9 
Assuming a standard deviation of .2 inches, has the average changed? 
SOLUTION: Hj: [f= 30.1 , Hy: fl¥# 30.1 , A= .05 
(1) Press [B'], 30.1, [R/s], .05, [R/s], 0, [R/s]. 
(2) Enter data using oa store () = .2 in R 


(3) Press 1 is displayed (reject Ho)» 
me ode) 0255472 47 ., 


07 


A new fad diet was tried out on 15 subjects and the Weight losses 
after one week were: 
Beeige 2.04 1.97 1.85 1.84 2.23 2.15 1.89 
1.93 1.99 1.86 1.90 2.09 2.16 2.04 
An advertisement claims that the weight loss after one week on the 
diet is at least two pounds. Do the data support the claim?(C¥= .05) 

SOLUTION: Hy? [kL m2 ; H, ii< Zz 

(1) Press [z'], 2, [R/s], .05, [R/s], +, [R/s]- 

(2) Enter data using iat 

(3) Press 0 is displayed (accept H), 

sl = 1302539836 . 

It is claimed that a certain drug will lower temperature within ten 
minutes. Five subjects having normal temperatures of 98.6 were given 
the drug and ten minutes later their temperatures were recorded in 
summary form as follows: }x = 491 , x = 48219 . Test the claim 


Biethe C¢ = .05 level. 


rs 








lL, 


5. 


SOLUTION: Ho: LL = 98.6 : H,: [L > 98.6 


(1) Press _eleiaier [R/s], 05, [R/s], ie [r/s] . 


as Cree 
(2) Store n= 5 in Ro3) x = 491/5 in R 9 and s =eEL= Ex in R 


1 
(3) Press O is displayed (accept HH)» 
ene me OOOO. 


x 14 


A student claims he can always answer more than half of the items on 
a true-false exam correctly, regardless of the topic. You devise a 
20 question exam on a subject of which he knows nothing and he 
answers 12 of the test items correctly. Would you conclude he has 
extraordinary powers? (C= .05) 

SCLUTZON: HH: eis aes Hy ' ee a 65 


Oo 
() ress [EF]. 5 R], 05, GS). BA. 


i)? Store nx = 12 in Ro, and n= 20 in Ro3' 


(3) Press O is displayed (accept Ho)» 
sl = .2517223358 . 
Of 694 respondents, 369 were in favor of more dissemination of birth 
control information. Is it safe to conclude that more than half of 
the population agrees with this position? (C¥= .05) 
SOLUTION: H Li a ee Ht i > .5 


(1) Press , 5, [B/s], .05, [R/s], 1, [R/s 


(2) Store nx = 369 in Ro, and n = 694 in Ro3 
(3) Press 1 is displayed (reject Ho)» 


xZt{ sl = .0474381802 . 
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L 
Oe 


Si. 


A soup can filling machine is supposed to fill each can with ten 
ounces of clear broth, with a variance of .01. A change in 
Waatability in either direction is undesirable. A random sample of 
20 cans yeilded )x* = 2107.7 and )x = 205.3 . Is the machine 
working within limits? ((¥= .05) 

SOLUTION: H:@*= .01 , : o* # .01 

(1) Press | 5' | ecole, [r/s], 0 5, |R/s|, 9, [r/s] . 

ire score n= 20 in Ro) yx =O a Roy? Ee = pCO main 


R and es Ozin 


O02 10 

(3) Press 1 is displayed (reject ny : 
A nail machine is supposed to manfacture i-inch nails with a 
standard deviation of .025 inches. A random sample of 30 nails 


yeilded a-esample value for s,, of .03 inches. Does this apparent 


A 


increase warrant shutting the machine down? ((C¥= .05) 
SOLUTION: Hj: * (025)? , H,: 02 = (.025)? 


(1) Press , .000625, [R/s], .0 eae Sleee [R/s] . 


(2) Store n = 30 in Ry, and (n-1)s* = 29(.03)* in R 


(3) Press QO is displayed (accept Ho) 


ieee 


A certain type of expensive electrical gear is supposed to have a 


1 


mean life of 1000 hours. The manfacturer is concerned if the mean 
departs in either direction from 1000. Five components were tested 
and they had the following burnout times; 

1075 1085 1060 998 995 


Is the mean still 1000? (CY= .05) 


fas 





SOLUTION: Hj: {£= 1000 , H,: ft # 1000 

(1) Press , 1600. ES, so5y [r/s], oF [R/s] . 

(2) Enter data using [ | 

(3) Press O is displayed (accept H_) 
Over a period of years there had been an average of 14 accidents per 
year in a certain city. This year the monthly totals were as 
follows: i 0 Ca tirce I at 3 ae il O 1 Cows 
Does this data agree with the theory that the number of accidents 
per month follows a poisson distribution with = t e67e (Cr= 205) 

SOLUTION: H+: \= 1.1667 , Hy: \# 1.1667 

(1) Press [z'l] , 1.1667, [r/s], 05, [R/s], 0, [z/s] . 

(2) Enter data using [> | 

(3) Press QO is displayed (accept HA) 

sl = .5703943172 
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PROGRAM 4 USER GUIDE - Two-Population Hypothesis Tests 
INTRODUCTION: The purpose of this program is to perform one-tailed and 


two-tailed hypothesis tests for the following two-population situations. 


BERNOULLI Py=Py for large sample sizes 

NORMAL Ly Ly for X and Y paired 

NORMAL [LeLly for X and Y independent with g,= gy*=c% 
NORMAL [i= [ly for X and ¥ independent with G)*, Gj, known 
NORMAL Uy= [Ly for X and Y independent with Cra 
NORMAL O;{=0;,," for X and ¥ independent 

NORMAL p = 0 for X and Y paired 


EXPONENTIAL NEY far X and Y independent 
GENERAL PROCEDURES: 

1. This program requires the Applied Statistics Module. After 
reading all three card edges, press to repartition (719.29). 

2. For each new problem press to clear all registers and to 
prepare for the following parameter entries. Enter C¥ and press 
then enter either 0, -1 or +1 for the desired alternate hypothesis, as 
in the table below, and press |R/S| . 

0 for Hy Bae oy Reieetaaledmtest 
-1 for Hy GC, < Oy , LoWer-tailed test 
+1 for Hy Oy >= Oy , Upper-tailed test 

3. For independent data entry (DEI Sequence) press | D |followed 
by data point Xs 1 Xa [R/s]. etc. for each x, Gai 2a) 
When all x, 'S have been entered press[ D jagain, followed by data point 


Ys? [z/s}, You? [r/s], etc. for each ve (GSA 2 aap lt) When all of the 


y;,'s have been entered press Ea! Mistakes in data entry should be 
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corrected immediately by reentering the unwanted data point and pressing 
[= +] , then enter the correct data point and press 
[z/s|. Alternate data entry is detailed in each subroutine. 

4, For paired data entry (DEP Sequence) press ioe followed by Xe 
[cet |, 5 [r/s], Xe [xzt],, ena [z/s], etc. when all of the paired 
data points have been entered press . Mistakes and alternate 
entry are as above. 

5. Each routine displays either a 1 for Reject ora O for 
Accept. The significance level is always in the T-register (R,) + 
PROGRAM 4 SPECIFIC PROCEDURES: 


Tests for BERNOULLI p,=p, for large sample sizes 
1. Press , enter CY, press [z/sI , enter 0, -1 or +i , and 
press [R/S]. 


2. Store n in Ryo Fee Rows m in B03 and my in Roy : 
9. Fress either 1 (reject) or 0O (accept) is displayed, then 


press to display significance level. 


Tests for NORMAL LL [Ly for X and Y paired 
Pee eress , enter CY, press , enter 0, -1 or +1 and then 
press R/S]. 
2. Enter data using DEP Sequence (Alternate entry: store (n-1) in 
Ro. and d(4/n/s,) in Ryo) 
Press either 1 (reject) or 0 (accept) is displayed, then 


press to display significance level. 


Ww 
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Tests for NORMAL [Ly = LL, with oe o ae =¢2 


i. 


Press {| B'{ , enter CY, press , enter 0 , -1 or +1 and then 

press (R/S) . 

Enter data using DEI Sequence (Alternate entry: store n in R 
: me =e 28 oe 

m in Ro3) x in Ry, y in Rog, sy~ in Rog and sy in Roo) 

Press either 1 (reject) or 0 (accept) is displayed, then 


press TO display significance level. 


iL 5) 


Tests for NORMAL fly = [ly OAg ae and ia known 


i 


3. 


Press , enter &, press RAF enter 0, -1 or +1 and then 
press ; 


Enter data using DEI Sequence and store Gre in R and Gis rl 


09 
Roo (Alternate entry: store n in Ryo) m in Ro3) x in Ry,» y in 


Ae Ce 
Rog? Oy im Rog and GY~ in Ro) 
Press either 1 (reject) or 0O (accept) is displayed, then 


press to display significance level. 


Tests for NORMAL UL, = [L, for Ce Z ge 


1. Press , enter (¥, press [R/s] , enter 0, -1 or +1 and then 


Za 


press : 


Enter data using DEI Sequence (Alternate entry: store n in Ry os 
. mae -. Din. 
min Ro3 eh Ry esi Rog: Sy in Rog 07): 


Press , i. (reject) or 0O (accept) is displayed, then 


and sy da 


press to display signiricance level. 
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Tests for NORMAL gj,” = Gy? 


Press , enter CY, press , enter 0, -i or +i and then 


De 


press ; 


Enter data using DEI Sequence (Alternate entry: store n in R 


: Le an 
m in Ro3) Ss, in Rog and s,~ in Ron): 


i: 


Press either 1 (reject) or O (accept) is displayed, then 


press to display significance level. 


Tests for NORMAL ae 0 


Le 


Press | E'! , enter CY, press Ee enter 0, -1 or +1 and then 


press : 


Enter data using DEP Sequence (no alternate entry). 


Press [ a’ | either 1 (reject) or O (accept) is displayed, then 


press to display significance level. 


Tests for EXPONENTIAL Nx = Ne 


Press , enter CY, press [r/s], enter 0, -1 or +i and then 


i. 


press : 


Enter data using DEI Sequence (Alternate entry; store n in Ryo) 


m in Ro3? x in Ray and y in Rog): 


Press either 1 (reject) or O (accept) is displayed, 


press to display significance level. 


PROGRAM 4 LABELS USED 


A 


A' SIN RST Ene 
B cos CLR xe 
e DEG RCL x 
D! INV STO - 
B GRD EXC 


8Q 


then 





PROGRAM 4 STORAGE REGISTER CONTENTS: 


2 2 = 
OO O, -1 or +1 07 Sy or Oy 144 x 
O1 my or o8 y er 
eZ az = 
wae yy“ OOM: oxy 16 - 24 used 
Seem or n 10 d(./n/s,) 25 (n-1) 
O4 nx or ie fal Gx 26 - 29 used 
O5 x" 12 used 
06 ) x¥ 13 used 


PROGRAM 4 SAMPLE PROBLEMS: 


1. 


According to last year,s statistics, 70% of the population were in 
favor of stricter smog control laws. This year only 65% of 1000 
respondents favored stricter laws. Does this represent a significant 
decrease? (CY¥= .01) 

SOLUTION: Ho! Py = Py ; Hy Py => Py 


(1) Press ee Oe [r/s] , +1, [z/s] - 


_ e a _ e 1 = 4 ° 
(2) Store n = 1000 in Ryo, nx = 700 in Ry, m= 1000 in Ro3 


and my = 650 in Rane 
(3) Press 1 is displayed (reject H,) 
sl = .0084920904 
Ten plots are split and half planted with variety A and half with 
variety B. The yields are shown below. Is there a difference in mean 
yield between the two varieties? ((C¥= .10) 
As 49 58 53 #60 45 %449 #66 55 AF 52 


Pees; 57 49 570 (Ue 853 
SOLUTION: oH: LL, ={f,» H, UL, # Ls 
(1) Press [s' | Por BAE 0, [R/s] - 


on 





B. 


(2) Enter data using DEP Sequence. 
(3) Press 1 is displayed (reject H,) 
sl = .0137674905 

A cigarette manufacturer tests tobacco from two different brands for 
nicotine content and obtains the following (in milligrams). 
Ai Bae 260. 25 ~~ 22 Zo 
B: ee? 28 Zo ee 26 
Do these results indicate there is a difference in mean nicotine 


content for the two brands? ((Y= .05) 


SMUONA MACE ft Hye he fl 

(1) Press ee OR [r/s]_, 0, [r/s] - 

(2) Enter data using DEI Sequence. 

(3) Press 1 is displayed (reject H.) 

esi. 01 70726612 

Two analysts make fifty independent determinations of the melting 
point of a certain chemical. The sample mean and variance of the 
data found by analyst I are respectively, 73.6 and 10 while the 
sample mean and variance found by analyst II are, respectively, 
72.4 and 8 . It is argued that there is a tendency for analyst I 


to get higher results. What is your conclusion? (C¥= .05) 


SOLUTION:  H: jL, = bly » Hy? [Ly > [ly 
(1) Press OR [r/s], ie [r/s] . 
(store mn — 50 in Rigs m= 50a Ro3) Sao Rays 


- 3 2a. 2 ey: 
y = 72.4 in Rog? Sy = 10 in Rog and Sy 8 in Rog: 


(3) Press 1 is displayed (xeject H,) 


sl = .022750062 


82 





We are interested in replacing wire B with wire A if the resistance 


per unit length is not significantly decreased. The results of 


twenty tests on each wire are presented below. If we know that the 


standard deviation of the two testing procedures are both .0017 ohms. 


What is your recommendation? (C¥= .01) 


A nO OL 
.049 
BO 
051 


O54 .049 .048 .049 049 .049 .053 .053 .049 
0047 048 .049 .049 +- «051 .051 1.050 .047 .050 


eeeeOohs.049 4051 OS 7 moot 6052 4050 1052 
oe O o>. .049  .052 OC OSODe, 852 .046 


SOLUTION: Ho! [Lb Lp , Hy Us > Lp ° 


(1) 
(2) 


(3) 


Press (coal sO, | 2/s|, LZ [ R/s |. 


Enter data using DEI Sequence and store Ors = (.0017)* 
2 oa Gk 
in Rog and Gy = (COL 77) ea Ron 
Press 0 is displayed (accept H, ) 
sl = .500837784 


Twenty plots of ground were planted with corn. Ten plots (Y) 


contained a special treatment. The variances were tested and found 


to be unequal. Is there a significant difference between the yields 


at the 5% level? 


2 = 13556 sy? = 02844 


‘ ‘ = , ok 
SOLUTION: Hj: [y= Ly» Hy: ly Ly - 


(1) 
(2) 


(3) 


Press [ st |, roy [2/s |, O, 


Store n= 80 an Ryo ne meer Ro3) = wom aN Ris ¥ = avis 


Ze ; a . 
x > 013556 in Rog and Sy 02844 in Roo: 


Press 1 is displayed (reject H_) 


Kocaeli = 0200 50012 


in Rog: § 


ry 


O 
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fee 


Suppose that two samples of ten and sixteen observations have, 


respectively, variances of .3888 and 2.25. At the 5% significance 


level Would you accept the hypothesis that Ge = on ? 
: : 2 Zz ; a Zz 


(1) Press OS. [r/s|, -1, [r/s|. 


(2) Store n= 10 in Rye 
2 


Sy 07° 
io) OQ is displayed (accept H,) 


x2t! sl = .994190704 . 


2 


meee 1th Ro3) Sy 


=e pooo an Rog and 


= 2.25 in R 


The following data measurements of students' ability in an IQ test 
are paired with scores in an acheivement test. Test for ei Gr, 
me «6105 954° 125 SE APO Oeil 90 132 116 
ee 47 4G 53 31 C4 43 75 40 80 5D 

SOLUTION: Hj: O=0, H,: PF 0 

(perress "| , .05, |R/s|. 0, |R/s\. 

(2) Enter data using DEP Sequence 

(3) Press 1 is displayed (reject H)) 

sl = .0012735515 

Two types of electric bulbs are observed as to length of life, with 
the following results: n = 46, m= 64, x = 1070 and y = 1041. are 


the mean lives significantly different? (CY= .10) 
SOLUTION! H+ Xx = Ay Ht AxyFry 


(1) Press [g'], .1, [R/s], 0, [R/s]. 


(2) Store n = 46 in Ry os cone = 070 4 


=) 
~ 
U 


14 and 


vy = 1041 in Rog : 
(3) QO is displayed (accept H,) 


sl = .5605970415 
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USER GUIDE FOR FIVE DISTRIBUTION PROGRAMS 

INTRODUCTION: The purpose of these five programs is to provide accurate 
approximate values for the following distributions: Normal, Binomial/ 
Multinomial, Chi-square, Student's t+ , and F., 

GENERAL PROCEDURES: Read in the appropriate program and proceed as 
below. 

NORMAL Distribution - Program 5 


Proceed with the following steps in any order. 


Enter Press Display 
1, Z A D(z) = normal density at z 
2 Z B Ea aje= 22 == 2) 
th : 
be «P(z) 6 ZA normal p percentile 
beg D a(z) = P(z =z) 
5. Z E NGA) Ee dr Fal 


BINOMIAL Distribution - Program 6 
Perform the first two steps and then do the remaining steps in 


any order. 


Enter Press Display 
1. n A n 
e P B Pp 
ae k é f(k;n,p) = binomial density at k 
oF kk D F(kjn,p) = P(Z=k) 
oe Ke 5 Q(k) = P(Z >k) 
6. - At pe 
7. - B' oe 
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MULTINOMIAL Density - Program 6 


Perform the following steps in order. Step 2 must be completed for 


each set (n,5P,) HOtMe, 2 cae, ee (k= 35) 
Enter Press Display 
Zi ns ae pay (can be used to check ae = N) 
Ps R/S ey (can be used to check Ps = 1) 


» - Eng otgy eo) 

* As an added feature the following capability exists in Program 6. 

4, store N in Rj, and press [B'] to display mw! (N S69). 
CHI-SQUARE Distribution - Program 7 

Store v in Ry and then the following steps in any order. 

Enter Press Display 
i a £(X°*) density at A 
2 ZZ 
2. x P(X?) = P(ZEX?) 
Z th , 
3. P(X *) xX p\¥) =p percentile 
* As an added feature the Gamma function can be evaluated as below. 


bey [D ] T (v/2) 


STUDENT'S t Distribution - Program 8 


Store v in Ry 5 and then the following steps in any order. 


Enter Press Display 


i |(Ct aqeeeedenaity ate 
a «(ot P(t) = P(Z=t) 
pa P(t) no = aes percentile 


* As an added feature the Gamma function can be evaluated as below. 


5. Store v in Ra 5 and press | p | to display [(v/2). 
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PeDistribution - Program 9 
Complete Steps 1 and 2 prior to Steps 3 or 4 btelow. 


Enter Press Dis ailtany, 


fle V x2 - 


oo 2 [A | ; 

3. F(¥,5¥5) P(F) = P(Z= F) 

4 , P(F) Be woe Es ae percentile 
PROGRAM 5 LABELS USED 

A B G D EB Dr RA 
PROGRAM 6 LABELS USED 

A B C D EB A' B' eu Be Oh 
PROGRAM 7 LABELS USED 

A B E D E A B' G D' E: 
PROGRAM 8 LABELS USED 

A B ¢ D EB B' G Ee 
PROGRAM 9 LABELS USED 

A B C D EB A' B' G D' Ey 
PROGRAM 5 REGISTERS USED 

feet 25 26 29 
PROGRAM 6 REGISTERS USED - ALL 
PROGRAM 7 REGISTERS USED 

fmeoge tO 1% 13° 44 415 417 thru 23 
PROGRAM 8 RHGISTERS USED 

Og 02 09 10 15 thr 24 
PROGRAM 9 REGISTERS USED 


00 02 08 09 11 15 thm 29 
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PROGRAM 5 - Sample Problems 
ier nid. (1.960). 


SOLUTION: 

Enter 1.960 and press| A | , 0 504409443 is displayed 
Mind P(1.960) = P(Z = 1.960) 

SOLUTION: 

Enter 1.960 and press| 3 _ |, 9750021748 is displayed 
Find P(-1.960). 

SOLUTION: 

Enter -1.960 and press| 3 |, 60249978252 is displayed. 
Find 4 ae 

SOLUTION : 

iter .05 and press| C |, -1.64521144 is displayed. 
Find 2 9s 

SOLUTION: 

Enter .95 and press[ c_|, 1.64521144 is displayed. 
Find Q(1.282) = P(Z>1.282) 

SOLUTION: 

Enter 1.282 and press| D |, .0999213886 is displayed. 
Find A(1.282) = P(zS |z| ) 

SOLUTION: 


Rover 1,252 and Press| E |, .8001 572228 is displayed. 
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PROGRAM 6 SAMPLE PORBLEMS 
Mmeerind £(k;n,p) given k= 6, n=10 andp= 4. 
SOLUTION: 
Enter 10 press then enter .4 press |B | then 
enter 6 press [ok .111476736 is displayed. 
2. Find P(8) when n = 10 and p= .75. 
SOLUTION: 
Enter 10 press then enter .75 press then 
enter 8 press on: 7559747696 is displayed. 
Peeeeind Q(2) when n = 10 and p= .25. 
SOLUTION: 
Enter 10 press then enter .25 press [3 | then 
enter 2 press PE |, 47407196 is displayed. 
4, Wind L and QJ for a Binomial distribution with n= 20 and p= .5. 


SOLUTION: 


Enter 20 press then enter .5 press then 
press ee 10 is displayed, then 
press (7* = 2.236067977 is displayed. 
oe Find F(n, ) for the following Multinomial case. 
ny 1 2 3 4 


Ds Pa el 


SOLUTION: 


Enter 10 press then 
enter n, P, [r/s| etc. and 


press , .00036288 is displayed. 
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PROGRAM 7 Sample Problems 


1. 


Find £(25) where v = 15. 
SOLUTION: 
Store 15 in Ras , then enter 25 and press 
.0134298528 is displayed. 
Find P(25) where v = 15. 
SOLUTION: 
Beore 15 in Ry 5 ,» then enter 25 and press 
9500565664 is displayed. 
: 2 
Find lie) ap 
in x 99 ) 
SOLUTION: 
Store 16 in Ry 5 , then enter .99 and press 
31 .99987803 is displayed. 
Find les). 
SOLUTION: 


Enter 16 press [ D |, 5040 is displayed. 


PROGRAM 8 Sample Problems 





ill 6 


Find £(2) where v = 20. 
SOLUTION: 
Store 20 in Ry 5 , then enter 2 and press 
.0580872152 is displayed. 
Find P(-.860) where v = 20. 
SOLUTION: 


Store 20 in Ry s , then enter -.860 and press 
.199990431 is displayed. 


a0 





3. Find + o9(14). 
SOLUTION: 
Store 14 in Ry , then enter .90 and press | 6 | 
1 .345266653 is displayed. 
4, Find | (4.5). 
SOLUTION: 
Store 9 in Ry 5 , then press oe 
11.6317284 is displayed. 


PROGRAM 9 - Sample Problems 





1. Find P(2.52) where v, = 5 and v, = 10. 


SOLUTION: 

Enter 5 and press rea then enter 10 and press 

enter 2.52 and press | 5], .8998498088 is displayed. 
2. Find P(.397) where v, = 10 andv , = 5. 

SOLUTION: 


Enter 10 and press lx2t |, then eter 5 and press 
enter .397 and press [ 3 |, .1002517964 is displayed. 
: ind 0,30). 
oe Fin F gae(t 0) 
SOLUTION : 


Enter 10 and press [eae then enter 30 and press 


enter .975 and press ie |, 2 oe 47 ts dvaplayed . 
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APPENDIX B 


SELECTED CHI-SQUARE INVERSE CDF APPROXIMATIONS 


DEGREES Fret Oe ye eed Goren, fl 7 
OF TABLED APPROXIMATE ACTUAL APPROXIMATE ACTUAL 
FREEDOM VALUE* PERCENTILE PROBABILITY PRRCHNGLoS rROBABITULTY 

©0039 0499 
1026 sooo 
©3531 .0 503 
ene MO) 5 Cu 
1.1460 60 500 
3.9404 10586 
PROG ~0 500 
10-8566 sO 500 
18.4927 10500 

same as TYPE [I 
PGW) 9001 
4.6052 9000 
Gee 57 39000 
(ae »9000 
9.2366 39000 
15.9873 .9000 
GAOL e «9000 
28 .4120 9000 
40.2560 .9000 

same as TYPE I 
7 8815 099 50 
10.5966 99 50 
125g 7 99 50 
14.8602 99 50 
16.7497 <7 oy 
25.1884 09950 
32.001 5 99 50 
O97 10 99 5D 
Eonorer 99 50 


* fret. Ze 465] 





S 2 


same as TYPE I 





SELECTED STUDENT'S t INVERSE CDF APPROXIMATIONS 


DEGREES 
OF 
FREEDOM 


TABLED 
VALUE* 


ie 
1. 
1. 
1. 
1. 
i 
1. 
1. 
i 
We 
i 


Lun SO 
as me a & 


rp TYPE [-— 


APPROXIMATE 


PERCENTILE PROBABILITY 


Rm RR RR RRR Ip kb 


LW Ein 0 ON 


os 


ACTUAL 





APPROXIMATE ACTUAL 
PERCENTILE PROBABILITY 
3249 .6000 
»2869 . 5994 
1274+ 5996 
2699 Sieh 
26oy Sy) 
»2660 . 5999 
27s .6000 
2507 .6000 
eo .6000 
250 .6000 
25 .6000 
=2909 .6000 
3.0777 »9000 
1.8764 .8993 
1.6383 9001 
Loe 7 9002 
del aZ 9002 
epi2o 9001 
1.3408 .9000 
1.3255 .9000 
pe.> .9000 
ico on .9000 
1.2958 .9000 
1.2886 £9000 
66 65 a7 soo 
9.7472 9948 
Bocce 99 50 
4.5990 099 50 
lL, 0299 .99 50 
Gog aZ 17) 
2.9476 19950 
2.8460 9950 
Zap t 199 50 
2.7047 199 50 
2.6604 99 50 
Zot 5 199 50 


pe eee i 





SELECTED F INVERSE CDF APPROXIMATIONS 


DEGREES 
OF 
FREEDOM 


TABLED 
VALUE* 


.0062 
0050 
0043 
0044 
0040 
0039 


eo 
7 


r—TYPE I— 


APPROXIMATE 


PERCENTILE PROBABILITY 


* [Ref. 2: pp. 472-U85] 


9.4, 


ACTUAL 


aoe eS) 


APPROXIMATE 


ACTUAL 


PERCENTILE PROBABILITY 








SELECTED F INVERSE CDF APPROXIMATIONS 


DEGREES r—- TYPEI— r—- TYPE II— 
TABLED APPROXIMATE APPROXIMATE ACTUAL 
VALUE* PERCENTILE PROBABILITY PERCENTILE PROBABILITY 


* [Ref . Za) DD. 472-485 | 





995 


WSEAS 62 99 50 
i Sos 9925 
Bi. Aes 99-44 
i Ooe 99 4 

9.361 99 4 
5/255 /9951 

Zoot Sal 99 53 
270 56 9964 
15.008 99 51 

Dole 99 50 
4.228 99 50 
3.59 9950 

25615.2 99 51 
285.06 9965 
13.228 9951 

4.070 199 50 
3.006 »99 50 
Cros 99 50 

ZCI 5 ec 9951 
ZO aoe 9965 
12.747 aoe 

Breer oo 
2.628 99 50 
1.984 99 50 

Be! sae 99 50 
289 .43 .9966 
12 197! 

Doe oo 
2-300 99 50 
1.606 99 50 
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